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ABSTRACT 

We use numerical simulations to model the effect of seeing and extinction modulations on weak 
lensing surveys. We find that systematic fluctuations in the shear amplitude and source depth can 
give rise to changes in the -E-mode signal and to varying amplitudes of large scale B-modes. Exquisite 
control of such systematics will be required as we approach the era of precision cosmology with weak 
lensing. 

Subject headings: cosmology: Lensing — cosmology: large-scale structure 



1. INTRODUCTION 

Weak gravitational lensing of distant galaxies by fore- 
ground large scale structure has emerg ed as a pow- 
erful tool for modern cosmology (see Mcllier 1999; 
iBartelmann fc Schneider! feOOli iRefregieil l2003al for re- 
views), which has already pro yided constraints on cosmo- 
logica l parameters (see e.g. iHoekstra. Yee. fe Gladdersl 
l200l Ivan Waerbeke fc Mellierl 1200.1 for the cur- 
rent status) and been touted for its potent i al to 
constrain dark en ergy llBenabed fc Bernardeaul 1 2001 ; 
Hutererl 12003 llul 12001 lAbazaiian fc Dodelsonl p003 : 
Benabed fe van Waerbeke! 120031 Uain fc Tavlorl 2003 : 
Heavens! 120031 lRefregierri2003b[ ITakada fc .TaTnl I2OO.I : 
Bernstein fc JaiiJl2004t iTakada fc Whitdl2004l l 



The technique relies upon the measurement of the dis- 
tortion that lensing induces in the shapes of galaxy im- 
ages. The percent level distortion induced by large-scale 
structure is generally referred to as cosmic shear. As the 
power of cosmic shear surveys increases, requirements 
that systematic errors in the measurement of these im- 
ages be accurately accounted for has become increasingly 
stringent. Fortunately nature has provided us with a 
means to test for some of these systematic errors. Since 
lensing arises from a scalar gravitational potential, the 
shear pattern it generates has a particular form. For 
example, the shear pattern around an isolated, spheri- 
cally symmetric mass distribution is tangential. Since 
this pattern has even parity it is often referred to as a 
(positive) _E-mode. In the absence of lens-lens coupling 
or higher order effects the shear pattern induced by lens- 
ing is pure E. A 45° rotation of the shear, to produce 
the (parity-odd) -B-mode, will null such a signal. Thus, 
a simple diagnostic test for a wide range of systematic 
errors is the presence of a -B-mode in the lensing maps. 

In this paper, we model the effect of systematic er- 
rors in the seeing correction and extinction on simulated 
weak lensing maps. We show that fluctuations in ampli- 
tude and depth of the signal across a field can generate 
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B-mode signals on large angular scales. However these 
B-modes do not closely track the change in the -S-mode 
power, and thus cannot be used to correct for these sys- 
tematic errors. 

We begin with a brief description of these effects in 
||5J and then describe how we generate simulated weak 
lensing maps in We discuss how we model the effects 
of seeing and dust extinction in the source population in 
21 and present our results in ^before concluding in ||SJ 

2. SEEING AND EXTINCTION 

One systematic effect that has received particular at- 
tention is the correction for the point spread func- 
tion (PSF) (e.g . iKaiser elHd1IT995t IHoekstra let, afl|l99l ; 
i KuiikerJ 119991 iKaiserl 1200(1 iBernstein fc Jarvisl 12002 : 



Hirata & Scliak 2003: Hoekstra 2004) and its anisotropy. 
Here, we consider two related issues that might prove 
troublesome: the effects of fluctuations in seeing condi- 
tions and of galactic extinction. 

Seeing causes a degradation in the lensing signal am- 
plitude by circularizing the images of background galax- 
ies. Corrections for this "isotropic" P SF effect have been 
studied using sim ulated images (e.g. iBacon et alJl2001l: 
lErben et, alJl200ll) and by marginalizing over the uncer- 
tain s hear calibration us ing a model fit to the power spec- 
trum ( Ish ak et "aT] [2004) . In practice a cosmic shear sur- 
vey consists of many pointings of a telescope, and seeing 
can easily vary by a factor of two between the best and 
worst pointings; even chip to chip variations in the de- 
tector can yield corrective factors that differ by as much 
as 10%. Systematic errors in the lensing measurement 
may therefore be introduced on the scale of the telescope 
pointings or of individual chips in the detectors due to 
imperfections in this correction. 

A second effect of seeing fluctuations is to modulate the 
effective depth of the survey by down weighting smaller 
images which have become too circularized. This not 
only reduces the signal to noise, it also alters the source 
redshift distribution, and can therefore be expected to 
have an effect that res embles source redshift clustering 
(e.g. !Bernardeaulll998l) . although on a different angular 
scale. The measured signal will then probe structures at 
different depths in different regions of the sky. Variable 
galactic extinction in a magnitude limited survey will 
have a similar effect, although it should not be correlated 
with the lensing signal. 
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3. SIMULATIONS 

In order to assess the systematic errors introduced by 
uncorrected changes in the calibration or depth of the 
survey we make use of simulated weak lensing maps. 
These maps are generated by ray tracing through N-body 
simulations. We use the method s and models described 
in detail in iWhite fc Valel l)2f)03j) . so we provide only a 
brief summary here. 

Our calculation is done within the c ontext of a ACDM 
model (model 1 oflWhite fc Vale 2003) chosen to provide 
a good fit to recent CMB and large scale structure data. 
The weak lensing maps are made from an N-body simu- 
la tion using a mu l ti-pla ne ray tracing code, as described 
in iVale fc White! (|2003). The code computes the 2x2 
shear matrix A 

At, -2 J " d X fl(x)ViV^ (1) 

where A describes the distortion of an image due to lens- 
ing, <f> is the gravitational potential, x is the comoving 
distance, and g(x) is the lensing weight 

,x(x s - x) 



g(x) 



dXs p{Xs 



Xs 
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for sources with distribution p(xs) normalized to J dp 
1. The shear matrix A is decomposed as 



72 ~ w 



72 + w 
K - 7i 



(3) 



where the 7, are the shear components, k is the conver- 
gence, and uj is the rotation, which is generally small. 

We make maps of the shear and the convergence at a 
range of source redshifts from z ~ to 3 in steps of A X = 
50 ft, _1 Mpc. In each case, a 2048 2 grid of rays subtending 
a field of view of 3° is traced through the simulation. The 
two shear components and the convergence are output 
at each source plane, and down-sampled to 1024 2 pixels. 
The final map is a weighted sum of the contributions from 
each source plane; for a distribution dp/dz s the weight 
given to source plane j is 



dp 

dZe 



H(zj)A X 



(4) 



We use a source distributio n of the form 
ijBrainerd. Blandford fc Smai]|IT99^ 

dp 



dz. 



oc z s exp 



-(z s /zo) 
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(5) 



For this distribution (z) = r(|)zo ~ 1.5 zq. We shall use 
zq = 2/3 for our base model, and include fluctuations in 
zq where appropriate. 

4. MODELING SYSTEM ATICS 

We are interested in the effects of modulations in the 
calibration of the amplitude of the shear signal and in the 
source redshift distribution. The former can arise from 
variations in the seeing correction for different pointings 
or for different chips, while the latter may occur as the 
result of seeing and extinction variations in a magnitude 
limited survey, as described above. 

We test the first of these effects by generating two 
"modulation maps" of unit mean. The first is a smooth 
gradient running from 0.9 to 1.1 from one side of the map 
to the other, while the second is generated by dividing 
the map into 64 22' x 22' regions, each of which is ran- 
domly assigned a value of either 0.9 or 1.1; this second 
modulation map has sharp edges, as we expect will oc- 
cur for chip to chip and pointing variations in the seeing 
correction. While 10% fluctuations are toward the up- 
per end of the range we expect, it serves to illustrate the 
effects of interest here. (A 5% fluctuation gives rise to 
.B-modes a factor of 4 smaller and i?-modes a factor of 
2 smaller than a 10% fluctuation, as expected.) These 
maps are then used to scale the shear at each pixel. 

We use a similar tactic to model the source redshift 
variation, using a third "modulation map" , also of unit 
mean, but now with a Gaussian distribution; unlike the 
shear calibration, we do not expect the seeing and dust 
extinction effects to contain sharp edges. We vary the 
amplitude of the modulation and the coherence scale of 
the fluctuations as follows. We first generate a Gaussian 
of zero mean and unit variance, with each pixel indepen- 
dent. We then apply a boxcar smoothing to the map 
with a scale of 4.5, 9 or 18 arcminutes. The map is then 
scaled to have 5, 10 or 20% standard deviation and we 
add 1 to every pixel. The map is clipped so that every 
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Fig. 2. — Variance in the aperture mass resulting from 10% modulations of the shear amplitude. (Left) The B-mode power (solid), and 
the absolute value of the change in the E-mode power resulting from modulations using a pixelized grid (dotted) and a smooth gradient 
(dashed) as described in the text. In both cases the i?-mode power is reduced on the largest scales for this particular realization. (Right) 
The B-mode power resulting from the same modulations. 



pixel lies in the range [0.1, 1.9] to avoid generating neg- 
ative or extreme source redshifts. The resulting map is 
then used to modulate zq in Eq. JSJ) when generating the 
shear maps from the sum of the planes described above. 
The shear at each pixel is thus a different weighted sum 
of the planes. In each case we find that the effect on 
different realizations of the shear field is very similar, so 
we show results from one 3° x 3° field below. 

5. RESULTS 

To quantify the changes in the E- and i?-modes caused 
by these modulations we compute the aperture mass 
statistics M ap and M± from the shear maps. The former 
should be sensitive only to i?-modes, while the latter is 
sensitive onl y to -B-modes. We use t he I = 1 form of M ap 
described bv lSchneider et ail l)1998j) . 

M ap (x ;R) = J d 2 x lT {x + x Q ) G (^pj (6) 

where xo is the position on the sky, R is the angular scale, 
7 T = —(71 cos 2(f> + 72 sin 20) is the tangential shear and 

G(y) = -L y 2 (i-y 2 ) fory<i (7) 

ttR z 

is the radial kernel, which vanishes for y > 1. On the 
smallest scales (< 2') the convolution is not well approx- 
imated by the sum over pixels, but this is not an issue 
for the larger scales which will be of most interest to us. 
To compute M±, we interchange 71 — > 72 and 72 — > —71 
before computing the integral in Eq. 

An example of an M± map, excluding the regions 
within R of the map edges, is given in Figure ^ 
Note that even the unmodulated maps contain some 
-B-mode power, as expected from effects such as lens- 
lens coupling and violations of the Born approx imation 
<|.Tam. Relink Whit.fll5ofifllVa.le Wnitell2003D . How- 
ever, modulations of the amplitude of the shear signal 
and of the source redshift distribution both significantly 
enhance the _B-mode signal. While the _B-modes are 
concentrated in regions where the amplitude or depth 
changes abruptly, there is additional structure in the 



lensing signal which makes the pattern somewhat com- 
plex. 

In order to quantify these effects further, we compute 
the variance of the M ap and M± maps, again excluding 
the regions within R of a map edge. This variance probes 
a narrow range of wavemodes in the power spectrum, 
peaked at roughly 1/3 of the filter scale, R. 

We show the effect of modulations in the amplitude of 
the shear in Figure [21 The smooth, gradient modulation 
generates very small J3-modes and a change in the E- 
mode power on all scales. This is not unexpected: the 
initial -B-mode is very small. If the transformation from 
shear to convergence (i?-mode) was completely local then 
rescaling the shear would not generate any -B-mode. The 
figure of merit is thus how much the gradient changes 
across the "non-local" scale in converting shear to con- 
vergence. Whatever the reason, the -B-mode is much 
smaller in amplitude than the change in the _B-mode, 
also shown in Figure [21 

A sharp modulation in amplitude gives a larger B- 
mode signal, with a similar change in the -E-mode sig- 
nal (Figure [21 • Below 25', the -B-mode signal in this 
realization does not track the change in the _E-mode sig- 
nal, which is decreasing for angles large compared to the 
modulation. However the amplitudes are more compa- 
rable than above, becoming very similar at the largest 
scales probed. 

The inclusion of fluctuations in the source redshift dis- 
tribution increases the large angle .B-modes by roughly 
2 orders of magnitude, as can be seen in Figure [3 The 
amplitude of this effect on scales larger than the modu- 
lation scale is (almost) independent of the angular scale 
of the modulation, but is an increasing function of the 
modulation amplitude in all cases. In extreme cases, the 
effect can be as large as 1% of the -E-mode variance (i.e. 
10% of the amplitude). The change in the -E-mode vari- 
ance also increases with increasing modulation, but not 
as significantly as does the B-mode. 

6. CONCLUSION 
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Fig. 3. — Aperture mass variance from modulations in the mean source redshift, using a coherence scale in the modulation map of 5'. 
(Left) The unperturbed _B-mode power (solid) and source redshift modulations of 5 (dotted), 10 (long dashed) and 20% (short dashed). 
(Right) The S-mode power, with line styles as for the E-mode. 



We have used numerical simulations to model the ef- 
fect of seeing and extinction modulations on weak lens- 
ing surveys. We find that fluctuations in both the shear 
amplitude and the source z-distribution can give rise to 
changes in the E-mode power and to large scale B-modes, 
with sharp changes in amplitude or fluctuations in the 
z-distribution giving the larger .B-modes. Since the B- 
modes do not closely track the changes in .E-mode power 
they cannot be used to "correct" for the above effects. In 
the case of strong S-modc enhancement, however, they 
can be used as a monitor for the effect. 

As we move from first detections into scientific ex- 
ploitation of cosmic shear, effects such as this will need 
to be carefully controlled. Photometric redshift informa- 



tion offers a likely route to mitigating fluctuations in the 
survey depth while offering many scientific advantages. 
The stable and excellent observing conditions from space 
can be expected to largely eliminate effects from point- 
ing and seeing fluctuations. Regardless of the route, it 
is clear that systematic errors such as those we have dis- 
cussed here must be controlled if we are to realize the 
full power of upcoming weak lensing surveys, which will 
help usher in a new era in precision cosmology. 

The simulations used here were performed on the LBM- 
SP at the National Energy Research Scientific Comput- 
ing Center. This research was supported by the NSF and 
NASA. 
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